Mechanism of hopping transport in disordered Mott insulators 
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By using a combination of detailed experimental studies and simple theoretical arguments, we 
identify a novel mechanism characterizing the hopping transport in the Mott insulating phase of 
Ca2-a;Sr2,Ru04 near the metal-insulator transition. The hopping exponent a shows a systematic 
evolution from a value of a = 1/2 deeper in the insulator to the conventional Mott value a = 1/3 
closer to the transition. This behavior, which we argue to be a universal feature of disordered 
Mott systems close to the metal-insulator transition, is shown to reflect the gradual emergence of 
disorder-induced localized electronic states populating the Mott-Hubbard gap. 

PACS numbers: PACS numbers: 71.30.+h, 72.15.Rn, 75.50.Pp 
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The interplay of electron correlation and disorder near 
metal-insulator transitions represents one of the most 
fundamental, yet least understood problems in contem- 
porary condensed matter physics. These issues are of par- 
ticular importance in strongly correlated materials. For 
example, recent studies on manganites have revealed the 
emergence of coexisting clusters near the M-I transition, 
which plays a crucial role in the mechanism for colos- 
sal magnetoresistance 0. Similar inhomogeneities have 
been observed in underdoped high-Tc cuprates, where 
spin glass and stripe formation feature disorder effects 
near the M-I transition . Despite extensive experimen- 
tal efforts 1', '5, 'a'l , however, the effects of the intrinsic 
randomness on the doping/substitution in Mott transi- 
tion systems remain largely an open problem, especially 
regarding the process of closing the Mott-Hubbard gap. 

Among the Mott transition systems, the single layered 
ruthenate Ca2-2;Sr2;Ru04 provides a rare case of the M-I 
transition that is highly susceptible to chemicalpressure 
because of strong lattice-orbital coupling [j, |y] . While 
one end member, Sr2Ru04 is a spin-triplet superconduc- 
tor , complete substitution of isovalent Ca changes the 
superconductor into the Mott insulator CaaRuO^by in- 
troducing distortions to RuOg octahedra ,8^ .il ilfl llll|. 
The lattice deformation in both pure and Sr substituted 
Ca2Ru04 stabilizes the insulating state by inducing or- 
bital polarization (orbital order) that creates a half-filled 
configuration near ep |3) Q> Q> Qi In addi- 

tion, the results by both Raman and optical conductiv- 
ity measurements indicate that the superexchange cou- 
pling J, and the effective bandwidth (cx VJ) in the in- 
sulating region are relatively insensitive to Sr content 
[l^ ITtI . Instead, the Sr substitution mainly controls the 
orbital polarization; the strong orbital polarization be- 
comes gradually relaxed with substitution, but suddenly 
vanishes at the M-I transition to stabilize the metallic 
state milllllllllli. 

This system is particularly suitable for the study of dis- 
order effects in closing of the Mott-Hubbard gap because 



the isovalent Sr substitution provides statistical random- 
ness without altering the key parameters: the net carrier 
density and the effective bandwidth. Furthermore, the 
variation in other parameters that usually complicates 
the analysis (crystal structure, phonon dispersions, etc.) 
should be small since a moderate substitution of Sr (< 10 
%) is sufficient to suppress the Mott insulating state. 

This paper presents the study of the effects of intrinsic 
randomness due to Sr substitution on the M-I transition 
in Ca2-2:Srj,Ru04. In this system, a well-defined Mott 
insulator appears at low Ts by the first order M-I transi- 
tion in < 2^ < 0.2. In the insulating phase, a systematic 
change in the transport and thermodynamic properties 
is observed, which elucidates the evolution of strongly 
localized states in closing the Mott-Hubbard gap. Fur- 
thermore, these localized states most likely reflect the 
tendency to spontaneously form inhomogeneities due to 
nano-scale phase separation in the vicinity of the Mott 
transition, as reported in other oxides Our analy- 

sis reveals that this situation gives rise to a novel mecha- 
nism for hopping transport, unique to insulating phases 
of strongly correlated systems with moderate disorder. 

Experiment. Single crystals of Ca2_a;Sra;Ru04 were 
prepared using a floating zone method jlSij . The molten 
zone was well stirred during the growth to mix Ca/Sr ho- 
mogeneously. X-ray powder diffraction show single phase 
samples without chemical segregation. Thus, Ca/Sr mix- 
ture provides a statistical disorder that prohibits classical 
percolation effects. Resistivity was measured by a stan- 
dard four-probe dc technique. Magnetization was mea- 
sured with a comlercial SQUID magnetometer. Specific 
heat was measured by a thermal relaxation method. 

Figure lA presents the temperature dependence of the 
in-plane resistivity Pah{T) measured on cooling . While 
Ca2Ru04 [x = 0) is a Mott insulator at low T and ex- 
hibits a M-I transition above 300 K shght 
Sr substitution dramatically decreases the M-I transi- 
tion temperature Tm-i to below 300 K. The transition 
is clearly seen by an abrupt increase of pab by a factor 
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FIG. 1: (A) T dependence of the in-plane resistivity p^hiT) 
measured on cooling for Ca2-a;Sri,Ru04 . Vertical broken 
lines are guides to the eye. (B) Phase diagram for the region 
< a:: < 0.2 with three different phases: paramagnetic metal 
(blue), paramagnetic insulator (red), and antiferromagnetic 
insulator (green). Tm-i and Tap are the M-I and AF tran- 
sition temperatures determined on cooling. Thick and thin 
lines indicate the first and second order transitions. 




FIG. 2; T dependence of the in-plane susceptibility for 
Ca2_a;Sr^Ru04 with x = 0.06, 0.09, and 0.15. AU curves 
were measured in a field-cooled sequence. The solid and open 
symbols indicate the results measured on cooling and heat- 
ing, respectively. Inset: (I) thermal hystereses observed in 
the paramagnetic phase for x = 0.06 and 0.09. (II) Cp/T vs. 
for X = and 0.15. The solid lines represent linear fit. 



of more than 10^. This change at Tm-i always involves 
large thermal hysteresis, indicating that the M-I transi- 
tion is first order. This is due to the first-order structural 
transition that occurs simultaneously with the M-I tran- 
sition 1^ . Reflecting the quasi- two-dimensional electronic 
structure, the anisotropy pciT) / pa\j{T) in the insulating 
phase reaches the value of the order of 10^ — 10** (not 
shown). As illustrated in Fig. lA and IB, pab in the 
metallic phase and Tm-i, both rapidly decrease with x, 
and the system finally becomes fully metallic aX x — 0.2. 

In the insulating phase, an AF ordering appears as in 
Ca2Ru04. Figure 2 shows the T dependence of the in- 
plane component of the susceptibility x(^) nieasured in 
a field-cooled sequence. The distinct increase in x{T) 
at low T arises from canted AF, as clarified by neu- 
tron diffraction measurements P|. The AF transitions 
are different in nature in the following two regions; in 
< X < 0.1, the AF ordering occurs at T < Tm-i with- 
out any hystersis, while in 0.1 < a; < 0.15, it coincides 
with the M-I transition. In fact, as for x = 0.06 and 
0.09, the M-I transition is reflected in the susceptibility 
hysteresis in the paramagnetic state (see the inset I of 
Fig. 2). The decrease of the paramagnetic x at Tm-i is 
attributable to the disappearance of Pauli component in 
the insulating phase. However, at a; = 0.15, a large hys- 
teresis is observed at Taf, indicating that both M-I and 
AF transitions at x = 0.15 occurs concomitantly with the 
first order structural transition. 

Here, we note two important facts that confirm the 
Mott insulating ground state in < x < 0.2. First, the 
AF ordering occurs only in the insulating phase. Second, 
p{T) does not show any anomaly at Taf for x = 0.06 and 
0.09. These facts agree well with the significant feature of 
Mott insulators: the separation between charge and spin 



degrees of freedom due to a large charge (Mott-Hubbard) 
gap compared to low energy spin excitations. 

Observation of hopping transport. In order to eluci- 
date the electronic state in the insulating phase, we an- 
alyze Pah{T) shown in Fig. lA. The result of Ca2Ru04 
(x = 0) fits well to activation type insulating behavior 

Pab(T) = ^exp(SG/2A:BT) (1) 

below 250 K to the lowest T measured, giving Eq ~ 
4500 K, consistent with the gap observed by the optical 
conductivity measurements |14LIi^. Since Ca2Ru04 is a 
Mott insulator, Eq should give the gap size between the 
upper and lower Hubbard bands (UHB and LHB). 
Variable-range hopping (VRH) conduction 

Pab(T) = ^exp(To/T)" (2) 

with a = 1/2 also describes the T dependence over al- 
most the same T region in agreement with Ref. [loj . 
VRH is usually observed in systems with strongly local- 
ized states near ep. However, fitting of the resistivity in 
a limited region alone does not give conclusive evidence 
for the presence of localized electronic states. In the case 
of Ca2Ru04, we measured the specific heat Cp{T) and 
found the electronic specific heat coefficient 7 to be 0±0.5 
mJ/molK^, as shown in the inset II of Fig. 2. Therefore, 
it is most likely that there are no localized states near ep 
and that activation-type behavior is the consistent inter- 
pretation of the low T transport in Ca2Ru04. 

To obtain the systematic estimate of Eq, we also try 
to fit the p{T) curve for each x to Eq. However, 
the well fitted region becomes rapidly narrower with x, 
indicating that the Sr substituted region does not obey 
activated behavior. In contrast, as in Fig. 3A, the VRH 
widely describes the insulating behavior. Equation Q 
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FIG. 3; (A) Inp vs. l/T^/^ and l/T^/^* for the insulating 
phase of Ca2-i,Sra;Ru04. The results for x < 0.09 are plot- 
ted against the lower horizontal axis l/r^/^, while these for 
a; = 0.15 are plotted against both and the upper axis 

Straight lines indicate a fitting result to Eq. (B) 
Corresponding plot for theoretical predictions. The results 
for W = 0.02 and 0.08 eV are plotted against the lower axis 
l/T^/^ while these for W = 0.16 eV are plotted against both 
l/T^/^ and the upper axis l/T^/^. Inset: (I) To obtained by 
fitting of p(T) to Eq. (2) with a = 1/2 for Ca2-:cSr^Ru04 
(solid circle) as a function of x (lower axis), and for theory 
(open circle) as a function of W (upper axis). (II) Model band 
structure used in our calculation. As the disorder increases 
{W : 0.02 eV 0.16 eV), the band-tails gradually fill-in the 
Mott gap, producing a soft gap similar to that predicted by 
Efros and Shklovskii, but on a much larger energy scale. 



with a = 1/2 describes Pah{T) for x = 0.06 and 0.09 over 
a respectable temperature range: 236 K-150 K (a; = 0), 
155 K-46 K (0.06), 152 K-14.0 K (0.09), and 21 K-7.3 
K (0.15). It should be noted that Pab(r) obeys VRH 
with a = 1/2 over a decade of T for a; = 0.09 and over a 
comparable range for x = 0.06. The x dependence of the 
fitting parameter Tq is given in the inset I of Fig. 3B, 
showing exponential decrease with x. 

For the highest Sr concentration x = 0.15, Pah{T) does 
not follow Eq. (0) well with a = 1/2, but better fits are 
obtained with a smaller value of the exponent a = 1/3 
over a decade of T (5 K - 65 K), as shown in Fig. 3A. 
The fit gives Tq = 8.8 x 10^ K. This value of the hop- 
ping exponent a is consistent with two-dimensional (2D) 
transport, within the Mott VRH picture Because 
our system has essentially a quasi 2D electronic struc- 
ture with large anisotropy pd Pah ~ 10'^ in the insulating 
phase, 2D hopping seems a natural choice. 

Origin of hopping transport. Generally, in the pro- 
cess of carrier doping in conventional semiconductors 
such as Si and Ge, the evolution of the transport behav- 
ior from simple activation to Efros-Shklovskii (ES) VRH 



(Eq. © with a = 1/2) and finally to Mott VRH has 
been observed, as the localized states fill in the band gap 
HH. The Mott VRH, observed at a; = 0.15 for a decade 
of T, usually appears when the system has a localized 
band with a nearly constant density of states (DOS) v{e) 
around ep- The low temperature Cp/T for x — 0.15 is 
well described by 7 -|- (3T^ with the same Debye temper- 
ature 0D = 410 K as Ca2Ru04 (411 K) and Sr2Ru04 
(410 K) (Inset II of Fig. 2). The electronic contribution 
7 is about 1.5 ±0.5 mJ/molK^. The extra entropy release 
of 0.36 J/mol K due to the increase in Cp/T below 14 K 
is much smaller than i?ln3 = 9.13 J/mol K expected for 
S — 1 moment at each Ru site, and is attributable to the 
freezing of the uncorrelated spins induced by disorder. 
Therefore, the 7 value above should give an appropriate 
estimate of v{£y) at T > 14 K. For 2D Mott VRH, the 
localization length ^ can be estimated using the formula 
fcBTo — 8.6/i^(eF)C^(i, where d is interlayer spacing. Our 
results yield a reasonable value of ^ « 25 A, that is sev- 
eral times the Ru-Ru interatomic spacing (« 3.8 A). 

In contrast, Efros-Shklovskii VRH usually appears 
when the long-range Coulomb interaction between lo- 
calized electrons is important, forming the so-called 
Coulomb gap (CG) with v{e) cx |e - epl in 2D [Uj. This 
leads to VRH with a = 1/2 and /cbTq ~ c^/k^, where 
K is a dielectric constant. However, Tg observed in our 
system (Inset I of Fig. 3B) is too large for long range 
Coulomb interaction. In fact, assuming k « 100 as in 
Ca2Ru04 PI, we obtain ^ of the order of 0.03 Afor x 
— 0.06, which is not physically meaningful. Correspond- 
ingly, the energy scale of fitting regions (T < 150 K) is 
also much too high in comparison with the ordinary size 
of CG (10-100 K) for semiconductors ij. It is therefore 
unlikely that the emergence of the Coulomb gap is the 
physical origin of the apparent ES-like behavior. 

What could then lie at the origin of the observed VRH 
with a « 1/2? To answer this question, we note that the 
gap-type structure in the DOS is essential for the expo- 
nent a — 1/2; The standard Mott VRH picture assumes 
an energy-independent DOS, leading to a small hopping 
exponent a = 1/3 in 2D, while the formation of the CG 
gives a strongly energy-dependent DOS, which in turn 
produces a larger hopping exponent a = 1/2. However, 
we have seen that the CG is too small to account for 
our observations. Instead, given the large energy scale of 
both the fitting region and Tq, the energy gap has to be 
0.1 — 1 eV, which actually fits well in the range of the 
Mott gap. Since the insulating phase of our system is 
a Mott insulator, we propose that the disorder-modified 
Mott gap lies at the origin of the VRH with a k, 1/2, as 
a new mechanism distinct from the ES scenario. 

Materials near a first-order transition separating two 
competing ground states are generally fragile against 
phase separation; moderate randomness can create co- 
existing clusters of competing ordered states 0]. In our 
case, the first order M-I transition separates the Mott 
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insulating phase and the paramagnetic metalhc phase. 
The statistical distribution of Ca/Sr inevitably creates 
(although small) the locally Sr-rich region with a broader 
local bandwidth. This statistical distribution of the lo- 
cal bandwidth should create extended tails of both UHB 
and LHB at their edges, which are easily localized. As 
the system approaches the metallic phase, the band tails 
extend deeper into the Fermi level region by enhancing 
their DOS. Finally, for large enough a:(Sr), these tails are 
expected to overlap and produce strongly localized states 
near ep. As the intrinsic disorder steadily increases, we 
thus expect the behavior to gradually cross over from ac- 
tivated, to ES-like VRH, and finally to Mott-like VRH. 
This is exactly what is observed in our experiment. 

Theory. To put these physical ideas on a more rigor- 
ous basis, we now present a simple theoretical model for 
hopping transport in disordered Mott insulators. First, 
we generalize the arguments of Mott and Efros- 
Shklovskii [ilj to a situation where the DOS has an arbi- 
trary energy dependence. According to these approaches, 
the transition probability between two localized states is 

a = aoexp{-2R/^-eo{R)/kBT}, (3) 

where Uo is an "attempt frequency," R is the distance 
between the two sites, eo{R) is their energy difference, 
and ^ is the localization length. The further the sites 
are from each other, the smaller the wave-function over- 
lap between the two states. However, sites with small 
energy difference are typically more distant from each 
other. The lowest energy So{R) of the accessible state 
within a radius R from a given site is implicitly given (in 
2D) by the solution of 

1 = TTRi^H{eo{R)), (4) 

where H{e) — v{e')de' . At any given temperature, 
the most probable hopping distance R is determined by 
minimizing the exponent of Eq. Q , giving 

1 ^ {-nH{eo)f^ 

e nu{eo)kBT ' ^ ' 

The solution of these equations determines R{T) and 
thus eo{T), which should be substituted in Eq. ©, to 
obtain the desired T-dependent resistivity p(T) ~ \/a. 
For the special cases of constant or linear forms for v{£), 
these equations can be solved analytically, reproducing 
the respective Mott (a = 1/3) or Efros-Shklovskii (a — 
1/2) hopping laws. In our case, v{e) assumes a more com- 
plex form, and a numerical solution is necessary. Given 
our fortuitous situation where the Sr substitution does 
not make significant change in net carrier density nor 
effective bandwidth, we chose a simple model for v{£). 
Specifically, we take the UHB and LHB to have constant 
DOS and bandwidth of 1 eV, with a gap of 0.4 eV, and 



f = 28 A, which corresponds to our material. The effect 
of disorder is described by introducing gaussian band- 
tails of width W, which tend to gradually fill-in the gap, 
as shown in the inset II of Fig. 3B. We examine the evo- 
lution of p{T) and the fitting parameter Tq to VRH with 
the exponent a = 1/2, as a function of disorder W . The 
theoretical results are presented in Fig. 3B for p{T) in the 
same fashion as the experimental ones in Fig. 3A, and in 
the inset I of Fig. 3B for Tq{W). The theory is able to 
reproduce all the qualitative and even some quantitative 
aspect of the experimental data for both p{T) and Tg, 
giving strong support to the proposed physical picture. 

In summary, detailed analysis of the Mott insulating 
phase in Ca2^a;Sra;Ru04 reveals that the Sr substitution 
to Ca2Ru04 changes transport from the activation type 
to the variable-range-hopping. The moderate random- 
ness intrinsic to the Sr substitution leads to the forma- 
tion of strongly localized states around the closing of 
the Mott-Hubbard gap. As a result, hopping transport 
emerges of the form similar to that predicted by Efros 
and Shklovskii, but with a distinctly different physical 
origin, and on a distinctly different energy scale. 
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